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ABSTRACT 

We present a statistical analysis of halo emission for a sample of 1047 edge-on disk 
galaxies imaged in five bands by the Sloan Digital Sky Survey (SDSS). Stacking the 
homogeneously rescaled images of the galaxies, we can measure surface brightnesses 
as deep as /Xr ~ 31 mag arcsec"^. The results strongly support the almost ubiqui- 
tous presence of stellar halos around disk galaxies, whose spatial distribution is well 
described by a power-law p cx r~^, in a moderately flattened spheroid (c/a ~ 0.6). 
The colour estimates in g — r and r — i, although uncertain, give a clear indication 
for extremely red stellar populations, hinting at old ages and/or non-negligible metal 
enrichment. These results support the idea of halos being assembled via early merging 
of satellite galaxies. 
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1 INTRODUCTION 

The most commonly accepted paradigm for the formation of 
structure in the universe, ranging from galaxies to clusters 
and superclusters of galaxies, predicts that they are assem- 
bled via the hierarchical clustering of dark matter (DM) 
halos, in the framework of the so called ACDM cosmology. 
Although these models have been very successful in describ- 
ing large scale structure and much recent progress has been 
made in understanding the processes behind the assembly 
of individual gala x ies th rough num erical simulations (e.g. 
Navarro fc White! Il994l; ISommer- Larsen. Gelato fc VedeM 
l^OTEtf^^Ts teinmctzl l200Ci IScarmaptec^^Tissere 



semi-analytic 



Kajjifmjrm ^_Wh^ Guidgrdonl 
Baugh. Cole fc FrenkI Il996l : 
iggih and their 



modelling 



ISomerville fc Primadd 

^ combination (e.g. 

J<auffmann, Nusscr fc Stcinmct j |l997|; iKauffmann et alJ 
1999,: Benson et al..,.200a: .Springel et alJl200ll) . much work 
is still needed. 

The study of the stellar halos of disk galaxies, such as 
our own Milky Way (MW), can give a substantial con- 
tributio n in this field. From many studies (for a review 
see e.g. iMaiewskillTggll . it is known that the MW halo 
is populated by old, metal-poor stars, whose origin is 
debated. They may have formed in the early stages of the 
dissipational collapse of the gas in the proto-galactic DM 
halo. Alternatively they may have been accreted through 
the stripping of stars from satellite galaxies. Evidence that 
this process plays an important role in building up the 
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stellar population of the halo, h as com e from the Sagittarius 
dSph (ilbat a. Gilmore fc IrwinI Il994^ ■ a low-latitude halo 
stream found in the Sloan Digital Sky Sur vey (hereafter 
SDSS, lYannv et all 120031: Ilbata et al.ll2003l) . the ongoing 
disruption process of t he Palomar 5 globular cluster 
JOdenkirchen et aljl20o3) and the detection of streams in 
local high velocity stars (|Helmi et alj|T999il . How much of 
the halo light and mass must be ascribed to this kind of 
interactions and how common they were in the past is still 
unclear. 

In order to answer these questions in a more general context, 
the extension of the observations to a statistical sample 
of external galaxies is required. Unfortunately the surface 
brightness contributed by a stellar halo similar to the one 
in the MW is typically 7 to 10 magnitudes fainter than the 
central parts of the other galactic components (disk, bulge) 
and than the sky. At present day, observations of the halo 
are available only for a handful nearby galaxies and pro vide 
varying results. Recent work bv lFerguson et abl \2Q0'2 ) has 
demonstrated the presence of substructure in the moder- 
ately metal-enriched stellar halo of M31, that is likely to be 
the relic of the disruption of one or more small companion 
galaxy. The search for halos in more distant galaxies, for 
which stars belonging to different components cannot be 
resolved, has been focused on edge-on disks, because of 
the much lower contamination of the halo by the projected 
disk stars. Problems with flat fielding, scattered light from 
other background and foreground sources, and extended 
PSF wings make it extremely difficult to obtain reliable 
photometry down to 29-30 mag arcsec"'^, as required in 
order to characterise galactic halos. A fter the first detectio n 
of halo around NGC 5907 claimed bv lSackett et al.l Jl994D . 
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many controversial results have followed. Th e latest and 
deepe st observations of this galaxy (see e.g. IZheng et alJ 
Il993^ strongly support the extraplanar emission being 
instead a ring which results from tidal disruption of a 
satellite galaxy. In the most complete study so far, (but 
see Sec. |3 for a more exhaustive review of recent results 
in the literature) , t he deep observations of a sam ple of 47 
edge-on galaxies, bv lDalcanton fc BernsteinI ll2002h . provide 
evidence for the ubiquitous presence of red envelopes 
around disk galaxies, which the authors attribute to an old, 
moderately metal-enriched, thick disk structure. 

In this paper we present the results from the statis- 
tical study of halo emission from a sample of more than 
1000 edge-on disk galaxies imaged by the Sloan Digital 
Sky Survey (SPSS. lYork et alJl200ol) . The SDSS is imag- 
ing about a quarter of the sky in the u, g, r, i, and z 
bands, with 54 sec drift scan exposures a t the dedicated 
2.5 m Apache Point Observatory telesco pe (lFulaigite_et_alJ 
| l996l:lGunn et al.lll998l: iHogg et al.ll2'00ltlSmith et alJl2002l: 
iPier et alJl2003^ ■ reaching ~ 25 mag arcsec at S/N ~ 1 
for a single pixel. In order to reach a surface brightness as 
low as 29 mag arcsec"^, we adopt a stacking technique, in 
which we combine the images of all the galaxies. First, the 
images must be geometrically transformed in order to make 
the galaxies superposable. Other external sources must then 
be masked. Finally, the count statistics of each pixel is con- 
sidered and a suitable estimator is chosen to represent the 
distribution. In this way it is possible not only to increase 
the S/N hy a, factor ~ VlOOO, but also to remove statisti- 
cally the major sources of contamination for deep photome- 
try of individual objects, namely foreground and background 
sources, inhomogeneities in the flat field and light scattered 
inside the camera. 

The paper is organised as follows: the sample is described 
in Sec. 13 in Sec. we describe the image processing and 
the stacking procedure. The results are analysed in Sec. [l] 
paying particular attention to the noise properties of the 
resulting images and to the possible sources of error and 
of bias. A discussion, including a comprehensive comparison 
with previous results in the literature and the possible impli- 
cations for galaxy formation, is given in Sec. El A summary 
and the conclusions of this work are reported in Sec. |S| 



2 THE SAMPLE 

As of April 2002 the SDSS has covered ~ 2000 square 
degrees both in 5 -band imaging and spectroscopy. 
iBlanton et alJ (|2002h have selected a sample of galax- 
ies from the Main Galaxy Sample dStrauss et al.l l2002l) . 
also known as 'Large-Scale Structure Sample 10' (LSSIO), 
including all sources with average surface brightness within 
the Petrosian^ radius in r band fi,. < 17.77, that have been 

^ iPetroaiai] il97(J) defines the Petrosian radius as the radius at 
which the surface brightness equals a given fraction of the av- 
erage surface brightness inside that radius. The Petrosian flux 
(and hence the magnitude) is defined as the flux insi de a certain 
number of Petrosian radii. See lStoughton et al.l i2002^ for the full 
description of the procedure adopted in the SDSS data reduction 
and the actual parameters. 



succ essfully targeted by the spectroscopic observations 
(see iBlanton et alJ l2003l for d etails about the 'ti ling' 
algorithm). We refer the reader to IBlanton" ietiOlliflQi) for 
the details of the sample and the regions of the sky covered. 
For the purposes of this work we selected a subsample of 
edge-on galaxies from the LSSIO requiring the following 
conditions to be satisfied: 

• Petrosian magnitude {Pmag) successfully measured at 
least in the three most sensitive SDSS pass-bands, namely 
the g, r, and i band; 

• Pmagi ^ 17.5; 

• isophotal semi-major axis^ a > 10 arcsec in i band; 

• isophotal axis ratio b/a ^ 0.25 in g, r and i band. 

Images of the 1221 selected galaxies were inspected by eye 
in order to prune from the sample objects that are unsuit- 
able for stacking. First of all we discarded a few percent 
of the galaxies whose axis ratios had been clearly underes- 
timated due to some failure in the SDSS PHOTO reduc- 
tion pipeline. Galaxies showing evidence of interaction with 
nearby companions, warps or other irregularities were re- 
jected. The absence of nearby bright sources contaminating 
the background was required as well. The resulting sample 
is composed of 1047 galaxies, ranging from —22.5 to —16.0 i- 
band absolute Petrosian magnitude, in units of mag -1-5 log h. 
The redshift distribution of the sample peaks at z ~ 0.05 
with a standard deviation of 0.035. Typical physical dimen- 
sions range from ~ 1 to ~ 25 kpc (Petrosian radius in 
the i band), with a median value of 7.4 kpc . A sum- 
mary of the sample properties is reported in Fig. Eye 
inspection of the single galaxies confirms that our sample is 
dominated by late- type disk galaxies, with at most a few per 
cents of Sb-c or Sb. Most of the galaxies appear to be almost 
bulge-free and the the remainder have quite small 'classical' 
bulges. 



3 THE IMAGE STACKING 
3.1 The image pre-processing 

The raw SDSS imaging data is available as 2048x1489 pixel^ 
(13.5x9.8 arcmin^) bias subtracted an d flat-fielded 'cor- 
recte d frames'. Using SExtractor v2.2.2 iBertin fc Arnout3 
|199G) we accurately determined the centre and the position 
angle of the major axis for each galaxy. From the segmen- 
tation image output by SExtractor we obtain a mask for all 
the sources in each field detected with S/N^ 1. The back- 
ground level for each galaxy was determined as the mode 
of the pixel count distribution in an annnlus with inner ra- 
dius 3 X aiao,i and outer radius 6 x aiso,i (where aiso.i is the 
isophotal major semi-axis at 25.0 mag arcsec"^ in i-band). 
This background was then subtracted from each image. 



2 the SDSS reduction pipeline PHOTO jLupton et al.1 l200ll) 

provides an eUiptical fit of the 25.0 mag arcsec"'^ isophote 
parametrised by the semi-axes a and b 
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Table 1. Median scale- lengths (in arcseconds) and rescaling ref- 




fkpc h-'l 




Figure 1. Distributions of the sample galaxies in Petrosian ra- 
dius (panel a), absolute i band magnitude (panel b), and redshift 
(panel c). 



3.2 Photometric scale-lengths 

The main hypothesis underlying the stacking technique as 
a reliable statistical tool to investigate halo properties is an 
assumed self-similarity of the disk galaxies, or at least of 
their diffuse component. There must exist a 'characteristic' 
scale-length, such that after rescaling, all galaxies can be su- 
perposed. The validity of this assumption will be discussed 
in more detail in Sec. 14.21 

Since the surface brightness distribution of edge-on disks 
is affected strongly by dust extinction, we chose the i band 
(Aeff = 7480A) as the reference pass-band in order to charac- 
terise the scale-length of our sample. This choice represents 
the best trade-off between the need for a high sensitivity 
and the desire to limit the effects of dust"^. Many different 
scale-lengths can be defined from photometric analysis. 
We consider the following four scale-lengths: 

• the Petrosian radius rpetro, as obtained from the 
SDSS photometric d atabase (see IStoughton et al.l 12001 
lAbazaiian et al.ll2003ll 

• the effective or half light radius rso 

• the exponential scale-length rcxp 

• an isophotal radius obtained from the one-dimensional 
light profile along the major axis riD 



and riD have been determined using dedicated soft- 



ware developed in the IRAF environment. In order to eval- 
uate the two former parameters, we derive a surface bright- 
ness profile from circular aperture photometry extending out 
to 1.5 isophotal (/li — 25 mag arcsec"^) radii, rso is then ob- 



tained as the radius enclosing 50 per cent of the total flux, 
while the exponential scale radius rexp is derived from the 
least squares fit to the surface brightness profile between r^o 
and rgo (i.e. the radius enclosing 90 per cent of the flux). 
The basic motivation for analysing an one-dimensional 
brightness profile, derived by collapsing the image of the 
galaxy along the minor axis, is to have an estimate of the 
brightness which is as much as possible independent of the 
inclination. The surface brightness enhancement is, in fact, 
important and severely dependent on the inclination, when 
the edge-on condition is approached. Assuming low or neg- 
ligible extinction, the quantity C, = —2.5 log where 
Fluxix) is the flux enclosed within the projected distance x 
from the minor axis of the galaxy, is independent both of the 
inclination and of the distance of the galaxy. As reference 
level for determining rio we chose C, — 25.0 mag arcsec"^, 
which corresponds roughly to a surface brightness of /i = 
25.0 mag arcsec"'^ for a typical edge-on galeixy. 
All four scale- lengths are quite strongly correlated, as ex- 
pected for self-similar exponential profiles. The typical scat- 
ter around the 1:1 relation between pairs of corresponding 
scale-lengths is ~ 15 per cent. 



3.3 The image transformation and stacking 

For each of the four scale-lengths considered, the median 
value in the sample, rounded to an integer number of pix- 
els*, has been taken as reference for the rescaling, as given in 
Table [3 All galaxy images have been geometrically trans- 
formed using the drizzle re-sampling technique (as imple- 
mented in the IRAF geotran task) in order to correctly 
propagate the original noise properties. Each image has been 
translated to the galaxy centre, rotated according to the 
measured position angle of the major axis, and expanded 
or contracted according to the ratio between the reference 
scale length and the corresponding measured one. The same 
transformation has been applied to the corresponding mask, 
and a 'grow radius' of 10 pixels (in the original image size) 
around each masked pixel has been adopted in order to mask 
the extended PSF wings of the brightest sources. Since each 
image has a different photometric calibration and the mea- 
sured surface brightness of a galaxy is differently affected by 
Galactic attenuation, depending on its location in the sky, 
an intensity rescaling is needed as well. The transformation 
from the original pixel intensity i to the rescaled one i' is 
given by the following formula: 

■ / ■ 1 r\ — 0.4( Zr,~ Z„ rf-f — X Airmass — Reddeninq) 

I — I X W ^ P ' , 

where Zp and Zp^ref are the flux calibration zero points of 
the image and of reference, respectively, K x Airmass is the 



^ The z band is in fact the red-most one, but its sensitivity is 
about a factor 2-3 lower than the i band. 



The pixel scale of the SDSS CCDs is 0.396 arcsec pix 
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airmass correction, and Redd ening is the attenuation due 
to ou r Galaxy, as given by ISchleeel. Finkbeiner fc Davis! 
l)l998h . All the quantities are expressed in magnitudes. 

The combination of the transformed and rescaled im- 
ages has been performed using the imcombine IRAF 
procedure. First of all, the masked pixels are rejected. The 
median of the count distribution of the remaining pixels 
is then calculated in order to obtain the median image. 
An average image is calculated as well, after clipping the 
16 per cent percentile tails of the count distribution of the 
unmasked pixels. Adopting the standard approximation for 
the mode of a distribution 



mode = 3 X median — 2 x average, 



we calculate the mode image. The resulting 512 x 512 pixel^ 
images extend much beyond the detected emission from the 
galaxies and allow us to determine the properties of the 
background in great detail. A careful comparison between 
the three statistical combinations shows that, despite of the 
conservative masking method adopted, the diffuse luminos- 
ity of sources other than the considered galaxies can result 
in significant skewness in the pixel count distribution and 
in a systematic increase of the sky surface brightness up to 
31 mag arcsec"^ in i-band in the outermost regions of the 
average frame with respect to the innermost ones, in which 
contaminating sources are avoided by selection. This system- 
atic effect in the background is almost completely removed 
in the mode image. In the following, therefore, we will al- 
ways refer to that one as the resulting stacked image. 
For each stacked image the residual background level is de- 
termined as the count mode in an annulus of 130 pixel inner 
radius, 96 pixel thick, and subtracted. 

In Fig. |5| we present the images for the four most sensi- 
tive pass-bands [g, r, i, z), obtained by stacking the images 
rescaled according to the galaxy exponential scale length. 
For the u band the signal-to-noise is not sufficient to say 
anything about the presence and the characteristics of halo 
emission, therefore we will neglect this band in the follow- 
ing analysis and discussion. Rescaling according to the other 
scale-lengths does not change the resulting images signifi- 
cantly, as we will demonstrate more quantitatively later on, 
and therefore they are not shown here. Intensity levels in Fig. 
121 are coded in square-root transformation grey scale, suit- 
ably adjusted to show the maximum extension of the low 
surface brightness envelope. We superpose isophotal con- 
tours in two magnitude ranges: black contours, correspond- 
ing to the faintest isophotes, are obtained adopting a boxcar 
smoothing scale of 10 x 10 pixel, whereas the white, bright- 
est ones, have been calculated with a 2 x 2 pixel smoothing 
scale. The represented isophotal levels are given in the figure 
caption. The images and the superposed isophotes show very 
clearly the presence of a diffuse luminous envelope around 
the disk, in all four considered bands. This halo (whatever 
physical meaning we give to this word) is particularly ev- 
ident and round in the r and the i bands, whereas it is 
significantly flattened in the g and z bands (but the lower 
sensitivity in the latter band does not allow us to say much) . 



4 THE ANALYSIS 

In this section we present the results of an analysis of the 
stacked images in the g, r, i, and z pass-bands. We mainly 
concentrate on the stacking of the complete sample (1047 
galaxies), because of the resulting higher sensitivity, but we 
also consider the stacking of three subsample of bright, in- 
termediate and low luminosity flat galaxies (Sec. 14.51 . in 
order to better understand the possible dependence of the 
observed halo properties on total luminosity. After analysing 
the background properties of the stacked images in order to 
assess our detection limits (Sec. 14. it . the photometric prop- 
erties are investigated by means of radial sector-averaged 
surface brightness proflles (Sec. 14.21 , and compared to dif- 
ferent models of luminosity distributions, namely thin-f thick 
disks and disk-|-halo (Sec. l4.3ll . The average halo colours are 
presented in Sec. 14.41 Unless otherwise specifled, we will al- 
ways use the stacking of the images rescaled according to the 
exponential disk scale-length. We devote an entire paragraph 
in Sec. 14.21 to demonstrate that the adopted scale-length is 
not critical in determining the characteristics of the observed 
halos. 

4.1 The background noise properties 

Since our photometric measurements are performed by in- 
tegrating the flux over different image areas, as will be ex- 
plained in Sec. 14.21 it is crucial to determine the noise prop- 
erties on different scales. The presence of large-scale fluctu- 
ations is already obvious from a superficial inspection of the 
contour plots in Fig. |21 We restrict the analysis of the noise 
to the annulus (130 pixel inner radius and 96 pixel thickness) 
on which the sky level has been computed, and we consider 
the rms of the intensity after re-binning the image by differ- 
ent linear factors, ranging from 1 to 50. The dependence of 
the rms on the scale-length L is well described by a power 
law of the form 

rms(L) = rmsQ L~°' , (1) 

where rmso is the noise on one-pixel scale, and a ~ 0.7 — 
0.8, thus indicating that the noise has significant large scale 
structure in excess to what is expected in the case of pure 
Gaussian noise (for which a — 1). 

4.2 The surface brightness profiles 

The most obvious way to investigate the properties of the 
approximately round, low surface brightness structures 
emerging in Fig. |5| is to extract surface brightness (SB) 
profiles, averaging the flux in large wedges at different 
position angles. First of all, we divide the image in four 
circular sectors of 60° aperture, centred at 0°, 90°, 180° and 
270° position angles (PAs). Each of these sectors is in turn 
radially divided into a number of coronae, geometrically 
spaced such that the outer radius of the feth corona is given 
by rk = ro ■ 1.15*, with ro = 6 pixel — 0.6 rcxp- Finally, 
we estimate the mean SB as a function of the radius for 
0° (i.e. along the disk) and 90° PA (i.e. perpendicular 
to the disk) averaging the SB in pairs of corresponding 
coronae in the two symmetric sectors at 0° and 180° and 
at 90° and 270° PA, respectively. In each of the four 
graphs of Fig. |3| we show the SB profiles for the g, r, i, 
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Figure 2. The images resulting from the stacking of the galaxies rescaled according to their exponential scale length, in g, r, i, and 
z band. Intensity is coded in grey levels through a square root transformation, adopting arbitrary bias and contrast. The overplotted 
isophotal levels are as follows: 30.0, 29.0, 28.0, 26.0, 24.0, 22.0 for the g and the r band; 29.5, 28.5, 27.5, 25.5, 23.5, 21.5 for the i band; 
28.0, 27.0, 26.0, 24.0, 22.0 for the z band. See text for the details. The scales indicate the offset from the galaxy centre in pixel, 10 
pixel=rexp. 



and z images respectively, in linear scale (upper left panel) 
and logarithmic scale (upper right panel): open circles 
represent the 0° PA profile, fiUed triangles the 90° PA. The 
shaded areas represent the level of the rms background 
fluctuations at the scale length corresponding to the area 
of the coronae in which the SB is averaged, as calculated 
in the previous section. Thus the 1 — a detection limits for 
the SB profiles can be assessed as ~ 31 mag arcsec"'^ for 



the g, r and i band and as ~ 28 mag arcsec ^ for the z band. 

In order to evaluate the errors on the SB we must con- 
sider contributions from 1) the background fluctuations and 
2) the intrinsic scatter in the average signal from a galaxy. 
The first one is just provided by the background rms as eval- 
uated in the previous section using Eq. with L given by 
the square root of the area of the corona over which the 
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Figure 3. The SB profiles in g, r, i, and z band, in 60° wedges along the disk (PA=0°, open circles) and perpendicular to it (PA=90°, 
filled triangles). Distances from the centre in pixels (10 pixels corresponding to 1 Texp) are shown in linear scale (left panels) and 
logarithmic scale (right panels). The sky rms is represented by the shaded areas. Dashed lines represent the 'pure' disk profile at r > 20 
pix, obtained as explained in the text. The bottom panels of each graph show the b/a axial ratio of the isophotes as a function of the 
major axis a. 



signal is integrated. The estimation of the intrinsic scatter 
contribution is much more cumbersome, since the shape of 
the statistical distribution of the pixel counts in the coro- 
nae is not known a priori and the hght distribution of the 
galaxy itself is likely to give a major contribution to the rms 
of the count statistics. In order to remove it, for each pixel 
we consider the deviation from the average intensity of the 
corresponding pixels in the four image quadrants, which are 
symmetric with respect to the x and y axes, and express the 
intrinsic rms as: 



\ 



1 

3~/V 



iV 4 



(2) 



where lij is the pixel intensity, A'^ is the number of pixels per 
corona per quadrant, and is the average intensity over 

the j index. In the absence of any large scale background 
fluctuation the error on the average intensity is just obtained 
dividing the ||5J by -\/4JV. Adding the contribution of the 
large scale fluctuations, we can write the error as: 



{rms{L))'^ 



+ 



AN 



(3) 



=1 j=i 



where rms{L) is given by Eq^with L given by the square 
root of the area of the corona, and the factor 2 comes in from 
having averaged two coronae. Such an error estimate has 
been proved to be consistent with the scatter of the average 
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intensity in the corresponding coronae in the four quadrants. 

In the lower panels of each graph in Fig. 13 we plot the 
axial ratio b/a of the isophotes as a function of the radius 
(semi-major axis). This is calculated as follows. We assign 
a number of SB levels and determine the radius r at which 
such SB's are reached at different PAs 9, by interpolating 
the SB profiles extracted from 6 wedges per quadrant. 
Assuming r — r(6) is well represented by an ellipse with the 
major axis along the x direction, we derive the best fitting 
semi-axes a and 6 by means of a standard least squares 
algorithm. 

As it is apparent from the straight-line behaviour of 
the open circles in the linear-scale plots of Fig. |3 the profile 
of the wedges centred on the disk show an exponential 
decrease of the SB with radius, as it is typical for disks. 
Perpendicular to the disk, instead, the rapid initial decrease 
of the SB becomes shallower as we go further away from 
the centre of the galaxy. The trend is well approximated by 
a power-law with index ~ 2, as the straight-line behaviour 
of the triangles in the log-scale plots points out. This is 
in general true for the four bands analysed. Ifowever, the 
relative intensity at the same radial distance along the 
disk and perpendicular to it is different in different bands, 
reflecting a possible dependence of the flattening of the 
diffuse halo on the band and, in turn, this points to colours 
gradients in the halo itself. As already noted from the 
isophotal contours in Fig. |5| the halo is prominent in the r 
and the i band, reaching a SB comparable to that of the 
disk at r ~ 8 — 10 roxp and making the isophotal shape 
significantly rounder {b/a ~ 0.5 — 0.6). In the g band the 
halo has very little extension (the surface brightness drops 
below 31 mag arcsec"^ beyond ~ 4 Texp), thus leaving 
the isophotes extremely fiattened. For the z band the 
measurements are inconclusive: there is some hint of a 
shallowing of the slope of the profile, but the S/N drops 
below 1 already at ~ 3 roxp at the level of 28.5 mag arcsec"^. 

In Fig. 13 for r > 20 pix, we also show as dashed 
lines the 'pure disk' profiles obtained subtracting the 90° 
profile from the 0° profile, after stretching the former 
by a factor 0.6~^ to take the flattening of the halo into 
account (see Sect. 1131 for details on the 0.6 factor). There 
is evidence that the 'pure' disk slope steepens beyond 
4 — 5 Texp in the r and i bands, consistent with studies of 
individu al disks in the literature (see e.g. I van der Krui3 
L2OOI; Kregel. van der Kruit fc de GriisI l2002t and refer- 
ences therein), whilst the results are unclear in the g and z 
bands. 

In order to rule out scattered light and the extended 
wings of the point spread function (PSF) as major contrib- 
utors to the observed halos, we have performed the same 
stacking procedure on the images of stars taken from the 
same frames as the galaxies to generate effective PSF's 
in each band. In each frame we identify a star 1) whose 
central brightness differs by less than 2 mag arcsec"'^ from 
the central SB of the galaxy and 2) which is distant from 
other contaminating sources, as requested for the galaxy 
selection. Then, sources other than the star are masked, the 
frame is re-centred on the star and the same geometrical 




Figure 4. The PSF profiles in g (solid line), r (dotted line), 
i (short-dashed line), and z band (long-dashed line), expressed 
as difference of magnitude with respect to the central surface 
brightness. The lines represent the analytic fit (Gaussian core plus 
exponential wings) to the radial profile measured on the stacked 
image of stars, as explained in the text. 



transformations as applied to the galaxy are performed. The 
radial proflle of the PSF can be reproduced by a Gaussian 
core plus exponential wings over a large extension. Fig. 
2] shows the analytic fit to the measured PSF's in the 
four bands, expressed as the difference of magnitude with 
respect to the central surface brightness. The decline in 
the core is very sharp, 6 — 7 mag arcsec"^ within 5 pixel, 
and the exponential wings contribute less than 1/10* of the 
central surface brightness at radii larger than 20 pixels. We 
will analyse the effects of the PSF on the measured profiles 
in the next section, by convolving with different models for 
the light distributions. 

The characteristics of the observed halo profiles show 
little or no dependence on the scale-length adopted for 
rescaling the images of the galaxies. In Fig. |5]we plot the 
SB profile in the 90° PA wedge in the four bands, using 
different symbols for the different rescalings: open circles 
for roxp, filled triangles for rso, open squares for risoiD and 
filled circles for rpetro. The r coordinate of each point is ex- 
actly rescaled in order to match the sample median values 
of the four scale-lengths considered. As in Fig.|3 the shaded 
areas represent the sky rms. The error-bars are calculated 
as explained before. The agreement between the different 
rescalings is extremely good: almost all points are consis- 
tent within the error-bars, with a handful of exceptions for 
the risoiD rescaling. Even for these, the deviation is less than 
3 (J. Thus we conclude that there is no systematic depen- 
dence of the average halo properties on the profile shape of 
the disk as described by variations in relative scale-lengths. 
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Figure 5. The surface brightness profiles in g, r, i and z band, at 90° PA, for diiferent adopted rescalings. Open circles refers to images 
rescaled according to rexp, filled triangles to rso, open squares to risoiD and filled circles to rpetro- The sky rms is represented by the 
shaded areas. 



4.3 Modelling the halo 

In this section wo investigate the structural properties of 
the detected emission by means of simple models of the 
light distribution. In particular, we consider the possible 
contribution from a thick disk component, with an ex- 
ponential vertical light density distribution, and from a 
(moderately) flattened power-law halo. Making predictions 
from models of this kind is non-trivial and must take into 
account that the observed emission results from a double 
convolution of the 'true' average light density distribution 
of the galaxies with 1) the distribution of inclinations 
and 2) the effective PSF, as computed above. The models 
are calculated as follows. First we assume a particular 
3-dimensional distribution of light. We produce a set of 
1000 Monte Carlo realizations of the projected surface 



brightness, uniformly varying the inclination angle of the 
disk between and 15 degrees (roughly corresponding to 
a projected axial ratio between and 0.25 for an infinitely 
thin disk, as required by the sample selection criterion for 
the galaxies). The 1000 realizations are then averaged and 
convolved with the analytic PSF computed in the previous 
section, separately for each band. We do not expect to 
reproduce the stacked images near the nucleus nor at 
small distances from the disk plane, because we do not 
model dust extinction. Thus we compare our models to the 
observations by means of vertical-cut profiles, that allow 
us to exclude these 'forbidden' regions. In each cut, whose 
width is chosen to be proportional to the distance from the 
minor axis of the image, we average the flux coming from 
the four quadrants. 
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In the thin+thick disk model the 3-dimensional hght 
density of each individual galaxy is assumed to be given by: 



-thin + R 



^thin 
( 

2thick 



(4) 



where A'' is the normalisation factor, R is the flux ratio of 



the thick over the thin disk, 



is the exponential scale- 



length of the radial SB distribution, zthin and Zthick are the 
exponential vertical scale-lengths of the thin and thick disks 
respectively. We fix rcxp = 10 pixels, and zthin = 1 pixel. 
No disk truncation is adopted. The realizations cover a log- 
arithmic spaced grid of four by four values in the zthick — R 
parameter space (zthick ranging from lO"'^ to lO'' * pixels, 
R from 0.01 to 0.4). The normalisation A'^ is left as a free 
parameter and fitted by minimising the x^- We find that 
the thin-|-thick disk models fail to reproduce the observed 
SB distribution of the observed halos, since they predict a 
sharp exponential cut-off of the SB along the z coordinate, 
which is inconsistent with the observed power-law profile. 
In the disk-|-halo models, we assume a halo component with 
a generalised Hubble density distribution as introduced by 
[eIsou. Fall & Freeman (1987), modified to allow the iso- 
density surface to be oblate spheroids; 



-(€+i)/2 



(5) 



where c/a is the fiattening parameter, rc is the softening pa- 
rameter or core radius^ , and ^ is the power-law index of the 
projected SB. The disk component is modelled by a double 
exponential distribution in r and z, without any truncation: 



(6) 



with rcxp and zo representing the exponential scale lengths 



in r and z respectively. We fix rc = 1 pixel and 



10 



pixels, and realize a grid of models in the h—^ — {c/a) param- 
eter space, covering the following ranges: h = 0.50 — 3.00, 
c/a = 0.3 - 1.0, C = 1-50 - 4.00. The total (disk-l-halo) 
normalisation and the normalisation of the halo relative to 
the disk for each model on the grid have been fitted by 
minimising the x^- The best fitting models for the g, r, i 
and z band are represented in Fig. |S| For each band we 
plot the vertical SB profiles at four different distances Xc 
from the centre of the galaxy, as obtained in the vertical 
cuts described above, offset by 1 mag arcsec"'^ one from the 
other. Even if the a-re extremely high, thus demonstrat- 
ing that the adopted models cannot reproduce in the details 
the complexity of the galaxy structure, the general agree- 
ment with the measured points is satisfactory and we are 
successful in reproducing the trend of the profiles. Besides 
the total model profiles (heavy lines), we plot also the ex- 
ponential disk components alone (light lines): while domi- 
nating at small height, they give negligible contribution at 
z > 30 pixels (3 rexp). The best fitting model parameters 
are reported in the panels of Fig. |S] The disk scale-height is 
quite well constrained to /i ~ -^rexp. The power-law slope 
(5) is steeper in g band (2.50), and increasingly shallower 
in r and i (2.00) and z band (1.50), while the halo shapes 



get increasingly rounder from c/a = 0.50 in g, to 0.60 in r 
and i, and 0.70 in z band. This is in good agreement with 
the previous analysis on the images themselves and on the 
radial profiles. The surface brightness of the halo component 
at r = rcxp along the minor axis is 25.84, 25.82, 25.03 and 
25.26 mag arcsec"^ in the g, r, i and z band respectively. 
We can estimate the amount of halo light coming from out- 
side the 25 mag arcsec"^ isophote as 2 — 3 per cent of the 
total galaxy light. 

The models derived above allow us to quantify the pollution 
of the outer envelope by scattered light and PSF wings. We 
note that, at r = 20 pix, the PSF-convolved disk component 
contributes ~ 30 per cent of the total surface brightness^. 
This contribution decreases very rapidly at larger distances, 
and becomes negligible at r > 30 pix in all the bands. We 
conclude that scattered light from the disk contributes much 
less than 30 per cent of the total (disk-|-halo) measured light 
beyond 20 pixel. 



4.4 The halo colours 

Based on the SB profiles presented in Sec. 14.21 we derive the 
colour profiles in the two 60°-aperture sectors including the 
disk and perpendicular to it. We concentrate on the g — r 
and r — i colours, excluding colours involving the z band, 
because of its lower sensitivity. In Fig. |7|the dots represent 
the colours measured perpendicular to the disk, with the 
error-bars derived from the errors on the SB profiles; for 
comparison, with the dotted lines we plot the colour profiles 
for the disk. This has a blueing gradient toward the outer 
parts, which is particularly evident for g — r, but it is still 
clearly apparent in r — i. Disk colours {g — r = 0.7 — 0.5, 
r — i =~ 0.4) are consistent with typical star-forming galax- 
ies, once we allow for significant dust extinction, as expected 
for edge-on conditions. The increasing presence of dust near 
the galaxy centre and age and metallicity gradients in the 
disk can also easily explain the observed gradients. Beyond 
4 — 5 rcxp there is evidence for colours getting redder, pos- 
sibly indicating that the halo is overtaking the disk. This 
is consistent with the steepening of the profiles of the 'pure 
disk', shown as dashed lines in Fig.|3 

The g—r profile perpendicular to the disk is extremely noisy: 
no clear trend can be established and the measurements 
are definitely unreliable beyond 20 pixels. The halo in these 
bands is roughly as red as the inner parts of the disk. The 
r — i colour shows instead a clear reddening toward the out- 
skirts of the galaxy, reaching extremely red colours ~ 0.8, 
about 0.4 mag redder than the disk. We derive the best es- 
timate of the halo colours using the mean colour around 20 
pixels and consider the uncertainties given both by the er- 
ror bars and by the scatter of the points around an 'ideal', 
smooth profile: 
g-r = 0.65 ±0.1 
r-i = 0.60 ±0.1 

These results are not corrected for the effects of the PSF. 
Based on the analysis in the previous section, we estimate 
that the SB at r = 20 pixel is polluted by scattered light 
from the inner parts by up to ~ 30 per cent. However, this 



^ The core radius is introduced only for mathematical conve- 
nience, to avoid the central divergence of a pure power law. 



^ 20 per cent in the g, 33 per cent in the r and the i, and 40 per 
cent in the z band 




Figure 6. Vertical cut SB profiles with superposed the prediction from the power-law halo+disk model. Points are the SB measured in 
four vertical cuts at distance Xc from the centre (in pixel units, see legend). Profiles are offset 1 mag arcsec"^ from each other. Heavy 
lines represent the model profiles, light lines are the disk component of the model alone. 



contribution is very similar in all the bands (20 per cent in 
the g, 33 per cent in the r and the i bands) and can red- 
den the real halo colours by some hundredths up to 0.2 mag 
(worse case for g — r). It is worth stressing that this as- 
sessment has large uncertainties, due to the very simplistic 
nature of the models here adopted and to their extreme sen- 
sitivity to the fitting parameters (e.g., changing the vertical 
scale-length of the disk affects very slightly the goodness of 
the fit, but can significantly change the disk and scattered 
fight contributions to the colours at 20 pixels). 
The robustness and the significance of the estimates given 
above, will be discussed in detail in Sec|S] 



tained from the complete sample. We notice that the halo 
SB is correlated with the luminosity of the galaxy, with the 
less luminous ones having also fainter halos. A similar cor- 
relation with the total luminosity holds for the disk SB as 
well, indicating that the relative brightness of the halo with 
respect to the disk is roughly constant. This luminosity de- 
pendence is barely observable in the g band, but becomes 
increasingly evident at longer wavelengths, with an average 
offset of ~ 0.5 mag arcsec"'^ between the 'bright' and the 
'faint' profile in the i band. This also implies that the halos 
of the bright galaxies are redder than the faint ones by some 
0.1 mag in both g — r and r — i colours. 



4.5 Dependence on the total galaxy luminosity 

In order to understand the dependence of the halo features 
on the galaxy luminosity, we analyse the stacked images (ex- 
ponential scale-length rescaling) for three luminosity bins in 
the i band, namely the 'bright' (—22.5 < iabs + 51og/i ^ 
-20.2), the 'intermediate' (-20.2 < i^bs + 51og/i ^ -19.0) 
and the 'faint' one —19.0 < iabs + 51og/i ^ —16.0). The 
profiles at 0° and 90° PA (obtained as in Sec. 14.21 are plot- 
ted in Fig.|S] open triangles for the 'faint' bin, filled squares 
for the 'intermediate', and open circles for the 'bright'. The 
solid lines running through the points are the profiles as ob- 



5 DISCUSSION 

The analysis performed in the previous section provides 
strong evidences for the generalised presence of a dif- 
fuse, low-surface brightness stellar component around 
disk galaxies. By stacking more than 1000 images, we 
are able to extract reliable photometry at SB level as 
faint as ~ 31 mag arcsec"^ in the g, r and i band and 
~ 28 mag arcsec"^ in the z band. 

It is worth stressing that the statistical estimator adopted 
in this work to compute the combined images is extremely 
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Figure 8. The stacked SB profiles in g, r, i and z band, at 0° (upper curves) and 90° PA (lower curves), for three luminosity bins: 
open triangles represent the 'faint' bin, filled squares the 'intermediate' and open circles the 'bright' one. The solid line is the all-sample 
profile. The shaded area is the rms background fluctuation. 



robust. The evaluation of the mode for each pixel in the 
composite image as 3 x median — 2 x average, with the 
average calculated over the count distribution after reject- 
ing the 16 per cent percentile tails and the masked pixels, 
is effective in removing outliers and spurious contributions 
from other sources, and in correcting for the skewness of 
the distribution, while reducing the noise of the standard 
mode estimation. 

The isophotal contours superposed to the images in Fig. 
121 show the presence of luminous halos whose shape is 
clearly rounder than the highly flattened central disk. The 
shape of the effective PSF, obtained from the stacking of 
1000 stars selected from the same frames and with similar 
criteria as the galaxies, after applying the same geometrical 
transformations, is not consistent with scattered light and 
extended PSF wings giving a major contribution to the 



light detected in the halo. Comparisons with the light 
distribution expected from simple thin-|-thick disk models, 
in which we have accounted for the whole range of galaxy 
inclinations and for the effective PSF, demonstrate that 
no exponential vertical disk component can yield the 
observed power-law shape of the SB profiles at r > 4 rcxp 
in any of the four considered bands. Adding a moderately 
flattened (c/a ~ 0.6), power-law (J oc r~'^ or p oc r~^) halo 
component to a simple exponential disk, it is possible to 
reproduce the basic features in the observed profiles quite 
well. 

It could be argued, that the emission we detect in the 
outer regions is produced by s m all bu lges. However, 
iMacArthur. Courteau fc HoltzmarJ i2003l) have shown 
that the bulges of disk-dominated galaxies, like the ones 
in our sample, are in general well represented by ex- 
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Figure 7. The g — r and r — i colour profiles. The filled circles 
refer to the 60°-aperture sector perpendicular to the disk, the 
dotted lines represent the colours in the wedge along the disk. 



ponential laws with /ibuigo/ftdisk ~ 0.13 ± 0.06 (/ibuigo 
and /idisk being the exponential scale-length of the bulge 
and of the disk respectively). Assuming an upper limit of 
fir ~ 18 — 19 mag arcsec"^ for the central surface brightness 
of the bulge, this implies that bulges contribute no more 
than 30 mag arcsec"^ beyond 20 pixels. The light we have 
measured is thus not a simple outward extension of the 
small bulges often seen in late- type galaxies. 
We have shown in Fig. |K| that the results obtained are 
independent of the scale-length adopted for rescaling the 
galaxies: this is a strong indication for the absence of 
any significant dependence of halo characteristics on the 
detailed shape of the disk SB profile, even if it does not nec- 
essarily imply complete homology of the disks in our sample. 

These results are consistent with some known features 
of the M il ky W a y hal o. Con firming previous results 
by iHarrij dlOTfill . IZinnI Jl985l) derived a halo steUar 
density declining as r~^'^ from the distribution of the 
globular clusters. Similar r esults were o btained from the 
RR Lyrae distribution by ISahal ^ll985^ ■ p oc r"^'^, by 
Preston^ Shectman fc Beer j lll99lH . p oc r-3-2±° i, and by 
Ivezic et alj j200Cf). o oc r~^-'^±0-'^ , The spatial distribution 
of B l ue Horizontal Branch (BHB) s tars p'rcston et al. 
Il99ll: iKinman. Suntzeff fc Kraf3 Il994ll is in good good 
agreement with those estimates, yielding p oc r~^'^ at 
height z > 5 kpc ~ 2.5 rexp above the disk plane. The SB 
provided by the MW halo stars is expected to be roughly 
of the o rder of 30 mag arcsec~^ a t r ~ 8 kpc ~ 4 rcxp in 
V band jBinnev fc Merrifieldlll993 . par. 10.5), co nsistently 
with o ur meas urements in q band. T he analysis bv lHartwickl 
^983) and bv lPreston et alJ lll99iri . considering RR Lyrae 
and BHB stars, provides an estimate for moderate flattening 



of the spheroidal halo c/a ^ 0.6, which is also consistent 
with the findings in this work. 

Comparisons with analogous studies of external galaxies are 
in general difficult. Although its small distance makes M31 
the easiest target for observing of the halo population of a 
disk galaxy, the prominence of its bulge makes it hard to 
disentangle the density distribution of the halo so that it is 
perhaps more appropriate to talk about a generic spheroid. 
iPritchet fc van den Berghl ^94) measured the SB of M31 
spheroid to fiy ~ 30 mag arcsec"'^ and concluded that 
it can be modelled either by a de Vaucouleurs law or, in 
its outer parts, by a power-law p oc r~^, which is much 
steeper than what we find. H owever, the gl obular cluster 
distribution follows p oc r" (iRacind [l99lh and there is 
evidence for a shallower power-law index in the outer parts 
of the halo from more recent observations (Irwin, private 
communication) . 

Recent observations of the red giant stars of the nearby, 
late-type spiral M33 (Ferguson et al. 2003, in preparation), 
seem to exclude the presence of a spheroidal component 
around this galaxy, but its nearly face-on aspect makes it 
difficult to draw firm conclusions. 

Results from more distant galaxies are even more un- 
certain, because of the overwhelming difficulties in going 
deeper than 28 m ag arcsec"^. After the first claim by 
ISackett et al.l ill 994) of the detection of halo emission from 
NGC5907, many discrepant measurements have been made 
by differ ent groups in bands from the optical to th e NIR 
(see e.g. iBarnabv fc Thronso"!] Il994l: ILeaueux et al.lll996l: 
Rudy et alJ Il997l: iJames fc Casahl Il998l: iLeaueux et alJ 
M9ar The latest observations bv Izh^i i^ (1^^) 
with intermediate-band filters, and by f^s^^D 1120001) 
in optical and NIR, togethe r with the RGB star counts 
derived by IZepf et alJ J2000l) from NICMOS observations, 
seem to rule out the presence of a halo, favouring instead a 
luminous ring produced by the tidal disruption of a dwarf 
companion. 

However, despite the non-detection of a ny diffuse compo - 
nent in the Scd galaxy NGC 4244 by iFrv et al.l Jl999t) . 
probably because of the low sensitivity of their obser- 
vations [hr < 27.5 mag arcsec"^), many studies during 
the last years support the idea of a luminous envelope 
(thick disk or halo) surrounding many of disk galaxies. 
lMorriso n_e^L|__^^97|) de tected thick disk emission from 
NGC 891: kbe et alJ EmSt) have measured R and / hght 
excesses with respec t to an ex p onent ial disk model in the 
Scd galaxy IC 5249; IWu et al.l (|20o3) have observed NGC 
4565 at 66OOA obtaining good accuracy photometry as faint 
as p = 27.5 mag arcsec"'^ and found a halo component 
with power-law r~^'^ Similar results are found by 

iRauscher et"ai] (|l99i) in the NIR K band for ESQ 240-Gll 
(power-law halo p oc r~^'^). After observing a sample 
of 47 extremely flat galaxies in B, R and K3 do wn to 
extremely faint SB, iDalcanton fc BernsteinI (JioO^) have 
claimed the ubiquitous presence of red stellar envelopes 
around disk galaxies. It is not clear whether these envelopes 
extend to a spheroidal halo or are just thick disks, as the 
authors claim, because their detection limits are just at the 
surface brightness where we start to see our power-law halo 
component unambiguously 

In order to establish the origins of these stellar populations, 
very precise colour measurements are needed. Due to the 
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Figure 9. Comparison of the halo colours with other stellar sys- 
tems in the (r — i) vs [g — r) diagram. The open circle with 
error-bars is our measurement, filled circles are the MW globular 
clusters from Harris (1996), the filled square with error-bar is the 
averageirms of the SDSS ellipticals at 2 = 0.05 from Bernardi 
et al. (2003), and the triangle is the measure of NGC 5907 from 
Lequeux et al. (1998). The shaded areas represent the g — r range 
for the thick disks in Dalcanton & Bernstein (2002) and the in- 
terquartile range of the Virgo and Coma ellipticals from GOLD- 
MiNe (Gavazzi et al. 2003, back hatches). 



enormous observational difficulties, reliable colour measure- 
ments of the halos around disks are very scarce: excluding 
the MW, M31 and M33 for which colour-magnitude dia- 
grams of halo stars and globular clusters can be obtained, 
the only optical colours avail able to now are the ones de- 



rived by|Leciueu2Let_al 



iDalcanton fc Bernstein 



199^ ) for NGC 5907. Beside these, 
20021) have measured {B — R) and 



{R — Ks) for the thick disks in their sample. The colours 
we derived in Sec. 14.41 are, in fact, very uncertain, both 
because of the intrinsic error of the measurements and 
because of the practical impossibility to avoid contributions 
from the disk component. The analysis in Sec. 14.41 showed 
that scattered light from the disk component is likely to 
affect the colours by a few hundredths of a mag, but colour 
excesses up to 0.2 mag cannot be ruled out. If we trust 
these results {g — r — 0.65 ± 0.1, r — i — 0.60 ± 0.1), we 
find that halos are made of stars which are only marginally 
(2(j) consistent with old, moderately metal-poor stellar 
populations. In Fig. |5]we compare the stacked halo g — r, 
r — i colours with different data taken from the literature^. 
Our point (open circle with error bars) is almost inconsis- 
tent with the colours of the MW globular clusters (filled 
circles, from iHarrisi il996il . being much redder in r — i, 



^ Photome tric transformation s between different standards are 
taken from lSmith et al.l j20o3). 



and as red as the metal-rich tail of the globular cluster 
distribution in g — r. We have highlighted 47 Tucanae in the 
plot (starred dot), because its integrated colours are very 
similar to MSl's halo. Its <; — r is about the same as our 
measurements. We note that g — r is also consistent with 
the blue end of the elliptical sequence as derived from the 
SDSS at the median redshift of the sample z — 0.05 (filled 
square, iBernardi et alJ 12003^. and from the observations 
in Virgo and Co ma by Gavazzi et al. in B and V, as given 
m GOLDMiNe (|g avazzi et al.ll2003tl . whose interquartile 
range is represented by the back-shaded area. Reasonable 
agr eement is found with th e optica l colours B — R derived 
by IDalcanton fc BernsteinI ||20o3) for their thick disks 
(whose range is represented by the shaded area in Fig.lill. 
Reconciling the measured r — i colour with any known 
stellar population is almost impossible, even if we allow 
for an extreme 0.1 mag reddening caused by the PSF, as 
discussed in Sec. 14.41 It is however i nteresting to no t e that 
the colours derived for NGC 5907 bv lLeaueux et alJ 
are consistent with ours, but nonetheless troublesome. We 
exclude significant dust reddening for two reasons. First, 
the unusually red colour in our data is r — i, despite the 
fact that dust mostly affects optical/blue bands. Second, 
in the r — i profile there is quite strong evidence for a red 
gradient toward the outer parts, whereas it is known that 
the dust is concentrated in a thin layer in the disk. Thus 
we conclude that there is evidence for the halos being made 
of extremely red stellar populations. This is likely to be 
primarily due to an old age, but other effects, such as high 
metallicity or 'exotic' low-mass dominated IMF, would be 
required in order to explain the unusually high r — i. 

Our red colours seem to exclude the possibility that 
the majority of the halos we observe around disk galaxies 
are made of metal-poor stars, or that they result from 
the integrated fight of globular cluster populations. The 
correlation between the disk and halo SB suggests a link 
between the two components. On the other hand, the high 
latitude extension of the emission rules out 'disk heating' 
as an effective formation mechanism. An accretion (or 
'cannibalism') scenario, in which the halo is built up by 
capture and disruption of spheroidal satellites, presents 
many advantages in explaining the observations. In this sce- 
nario dwarf spheroidal galaxies, made of old, metal enriched 
stars, are tidally disrupted by the gravitational field of the 
central galaxy and their stars are spread to fill the phase 
space almost isotropically in a few dynamical times. Thus 
both the spheroidal shapes and the extremely red colours 
of the halos could be explained, at least qualitatively. As 
already mentioned, recent observations in the halo of MW 
and M31, of the Sagittarius stream in the MW, and of the 
ring in the halo of NGC 5907 support the idea that this 
mechanism has been working till very recent times and is 
not uncommon among the disk galaxies. 

Our analysis in different luminosity bins shows that 
the halo luminosity, at least on average, is proportional to 
disk luminosity. Prominent halos are in the more luminous 
galaxies. The decrease of the average surface brightness at 
lower luminosity affects mostly the redder pass-bands and 
is reflection of the well known correlation of disk surface 
brightness with disk luminosity in the red bands (see e.g. 
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IShen et alJl2003|) . The light we see does not appear to be a 
straightforward extension of the bulge, since the power-law 
shapes of the profiles, irrespective of luminosity, are not 
consistent with classical de Vaucouleurs or exponential 
bulges. We can interpret the halo-disk luminosity relation 
we find here in the hierarchical picture, in which more 
luminous galaxies sit in more massive DM halos, with 
a larger number of merging sub-halos contributing to 
the stellar halo luminosity. However, deeper and more 
detailed observations, along with more reliable theoretical 
predictions for the number, stellar mass and metal content 
of the accreted satellites, are needed in order to confirm 
this hypothesis. 



6 SUMMARY AND CONCLUSIONS 

By stacking a large number (> 1000) of edge-on disk 
galaxies imaged in the SDSS we have been able to detect 
a diffuse, spheroidal, low-surface brightness component 
around the disk. This detection is significant in the g, r, 
i and z bands, and cannot be ascribed to any obvious 
instrumental artifact (e.g. scattered light or PSF). Given 
the statistical estimator we adopt for combining the images, 
our result indicates that a substantial fraction of the stacked 
galaxies must share the observed halo characteristics, even 
if we cannot exclude the possibility that a number of disk 
galaxies actually have no halo at all. The halo can be 
described by a power-law projected density profile I oc 
with a ~ 2 nearly irrespective of the band. 
The colour measurements provide inconclusive and trouble- 
some results, but there is a clear indication for extremely 
red colours, g — r is consistent with old, moderately metal- 
poor stellar populations, such as the more metal-rich MW's 
globular clusters, 47 Tucanae, the halo of M31 and the most 
metal-poor ellipticals, r — i is (at 2 cr) 0.2 mag redder than 
the reddest known stellar populations in globular clusters 
and elliptical galaxies and it is difficult to reconcile with 
any theoretical models, even allowing for ad hoc modified 
IMF's dominated by low-mass stars and high metallicity. 
The data also suggest a correlation between the luminosity 
of the halo and the total luminosity of the galaxy. 
The results presented in this work are far from being 
conclusive, but nevertheless they are consistent with the 
idea that a large fraction of disk galaxies are surrounded 
by a luminous halo. The colours, although affected by large 
uncertainties, hint at old, but not particularly metal-poor 
stellar populations, thus supporting a scenario in which 
the halos are mostly contributed by stars stripped from 
accreted or merged companions, in which the chemical 
evolution was already advanced. 

Deeper, individual observations of a large sample of nearby 
galaxies will be required, however, in order to assess the 
validity of this scenario by quantifying not only the average 
halo, but also the whole distribution of individual halo 
parameters. 
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